Bulk elastic fingering instability in Hele-Shaw cells 
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We demonstrate experimentally the existence of a purely elastic fingering instability which arises 
when air penetrates into an elastomer confined in a Hele-Shaw cell. Fingers appear sequentially 
and propagate within the bulk of the material as soon as a critical strain, independent of the elastic 
modulus, is exceeded. Their width depends non-linearly on the distance between the confining glass 
plates. A key element in the driving force of the instability is the adhesion of layers of gels to the 
plates, which results in a considerable expense of elastic energy during the growth of the air bubble. 



PACS numbers: 61.41. +e,62.20.-x,68.35.Np, 68.35.Gy 

Bulk fingering instabilities in viscous liquids confined 
in Hele-Shaw cells, commonly known as the Saffman- 
Taylor instability, have given rise to considerable exper- 
imental [1, and theoretical effort [2^. In the context of 
the liquid to solid transition (gels [3], foams [4], yield 
stress fluids and Maxwell liquids [8^), this instabil- 
ity translates into a fingering to fracture transition. 

Purely elastic instabilities in a confined geometry have 
received much less attention. An elastic fingering insta- 
bility has been observed during the peeling of a thin 
layer of elastomer from a rigid substrate O [10] or in 
probe tak experiments [11] where a semi-spherical probe 
in contact with the soft solid is pulled up at a con- 
stant speed. In most cases, the instability settles on the 
crack line, which is observed to progress at the interface 
between the elastomer and the rigid body it separates 
from [9l [12] . The energetic balance involving adhesion, 
elastic energy stored within the film and rigidity of the 
plate favors a meandering instability of the crack front, 
which has been fully characterized experimentally and 
theoretically |13| [14j. In contrast, Shull et al. [fSl [T6] 
have observed fingers propagating within the bulk of an 
acrylic tri-block copolymer gel. They demonstrated the 
reversibility and the non hysteretic character of the pro- 
cess. However, to our knowledge, there has been no quan- 
titative characterization of the observed patterns yet, and 
one dramatically lacks experimental data to constrain the 
possible mechanisms at the origin of this instability. 

In this letter, we demonstrate the existence of a large 
scale - centimetric - elastic fingering instability arising 
within the bulk of a polyacrilamide gel confined in a Hele- 
Shaw cell. The instability appears when a critical strain, 
independent of the shear modulus is exceeded. Fingers 
appear sequentially, and, ultimately, lead to the spectac- 
ular fiower-shaped pattern displayed in Fig. [l] The size 
of the fingers increases non-linearly with the cell thick- 
ness, while the number of fingers decreases sub-linearly. 
The latter is independent of the shear modulus. 

Materials and methods. We use polyacrylamide gels 
made from acrylamide monomers and bis-acrylamide 




FIG. 1: Elastic fingering instability. Pattern of the desta- 
bilized interface {G — 650 Pa; b — 2.1mm; Q — 126ml/min) 

cross-linkers [17]. The relative concentrations of these 
two components control the shear modulus of the gel 
from G = 80 to 1060 Pa. The strain to fracture of the 
produced materials decreases with G, but exceeds 300% 
in all cases. Adding dye to the gel allows to see where 
the material sticks to the glass plates. Cyclic shear rhe- 
ology of the gels with dye indicates a very small ratio of 
the elastic modulus to the loss one: G" /G' < 10~^ for 
frequencies ranging from 0.01 Hz to 100 Hz. 

Two different experimental setups are used. Setup 1 
(Fig.[2]left-(a)) is a classical Hele-Shaw: it consists in two 
10mm thick glass plates of lateral sizes (250mm x 250mm) 
separated by thin spacers of thickness h G [0.5 — 5] mm. 
Setup 2 (Fig. [2]left-(b)) is an original design, which con- 
sists in a sealed Hele-Shaw cell made of two 10mm thick 
glass plates of lateral sizes (250mm x 125mm), with two 
opposite mobile walls, acting as pistons, and pulled at a 
prescribed velocity V by synchronized step motors. The 
gap is fixed to the value h = 2.1mm. This setup is closer 
to classical tensile tests in solid mechanics. In both ex- 
periments the cells are filled with polyacrylamide before 
gelation. During this process, we maintain an initial air 
bubble of controlled diameter (I^o = 23mm in Setup 1, 
1^0 = 15mm in Setup 2). In Setup 1, air is blown with 
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FIG. 2: Setups and basic observations: (color online) (Left): (a): Setup 1, classical Hele-Shaw cell with air injection at 
the center, (b): Setup 2, sealed Hele-Shaw cell with two moving walls acting as pistons. (Center): Sketch of the instability 
observed in Setup 2 {G — 500 Pa, 5 = 2.1mm): Evolution of the perimeter P of the pattern as a function of time, exhibiting 
four distinct stages: (1) Circular growth; (2) Nucleation of the fingers; (3) Swelling of the fingers; (4a) Crack appearing in the 
gel layer leading to (4b) an interfacial fracture. Pc (respectively Pr) is the bubble perimeter at the onset of fingering (reps, just 
before fracture). (Right): Reversibility and hysteresis observed in setup 1 (C = 840 Pa, 6 = 1.4mm): Pressure 11 vs. strain 
7 = (P — Pi) /Pi, (red curve : inflation; blue curve: deflation). 



a syringe pump in order to grow this bubble. The pres- 
sure varies from to 1.5 bar. In Setup 2, the depression 
induced by the motion of the pistons sucks air into the 
gel. Because polyacrylamide is incompressible, the speed 
and size of the pistons impose an air flow rate equal to: 
Q = 2WbV = 126ml/min. In both cases, the loading 
rates are such that the material can safely be considered 
as purely elastic. 

Most observations are top views of the cell recorded at 
50Hz for Setup 1 and 1500Hz for Setup 2, which gives 
an idea of the robustness of the instability over a wide 
range of interface speeds. We measure the thickness of 
the gel layers through light absorption. When perform- 
ing such measurements, we use a silicon oil (Rhodorsil 
V20), which is immiscible with the gel but has a similar 
refractive index and therefore avoids refraction. The ob- 
served patterns show no differences with the one obtained 
when injecting air: The interfacial tension between the 
injected fluid and the gel is not a relevant parameter. 

Experimental results The same scenario holds in both 
setups and can be summarized as follows (Fig. |2]-center 
and Supp. Materials [24 ). The evolution of the perime- 
ter P of the pattern allows to distinguish four stages. 
During stage (1), a circular bubble grows. Shades of 
gray behind the dark line of the front show the existence 
of layers of gel adhering to the glass plates. In stage 
(2), fingers burst out successively, with an experiment- 
dependent order of appearance, and with a speed much 
larger than the velocity of circular extension. This leads 
to a sharp increase of the perimeter. Stage (3) corre- 
sponds to the swelling of the fingers once they are all 
formed, and the increase of the perimeter is much slower 
again. Then, in stage (4) the layer of gel behind the tip 



breaks, leading to an interfacial crack. 

Prior to the interfacial fracture stage, the phenomenon 
is completely reversible: Fingers deflate when the ten- 
sile stress is released, and the initial bubble is recovered. 
The reversibility of the process clearly demonstrates the 
purely - possibly non linear - elastic nature of the insta- 
bility. Figure [2] (right) displays the pressure 11 as a func- 
tion of strain 7 = {P — Pi)/ Pi^ where Pi — ttDq is the 
initial bubble perimeter. It shows that the phenomenon 
is slightly hysteretic. When air is injected a second time, 
the pressure follows the red curve again, which discards 
the hypothesis that the material has undergone damage 
during the first cycle. 

By filming the cell at an angle, we could see a clear 
meniscus on all fingers. Besides, as already shown in 
Fig. [2| center), interfacial cracking occurs at a late stage, 
after the full development of the instability. There two 
distinct fronts are easily observable, one corresponding 
to the interfacial crack, the other one to a deformation of 
the bulk ahead of it (Fig. [2| center) image 4(b)). Further 
quantitative evidence of this 3D character is provided 
in Fig. [3j where color encodes the gel thickness mea- 
sured from light absorption (Fig. |3^-b). Before a finger 
bursts out, the interface is already locally deformed. The 
gel thickness increases with the distance to the injection 
point. The finger first develops with a constant width and 
almost with a constant thickness as further illustrated in 
Fig. [sj^c), where the finger profiles within the thickness 
of the cell have been plotted at constant time intervals, 
assuming a symmetric shape. The finger grows up to a 
certain size where only the extremity is swelling while 
the back, forming a saddle, does not evolve any more. 
This early stage is followed by a regime where both the 
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FIG. 3: Formation of a finger: (color online) (a):Top view 
before (left) and after (right) the finger has formed. The 
initial front (black dotted line) recedes in the finger vicinity. 
Color bar: thickness of the gel in b units, (b) : Thickness 
profiles of a finger developing in time. Color codes 137 time 
steps from blue to red separated by 0.6 ms. G = 100 Pa. 



extremities and the "saddle" parts of the fingers increase 
in size. Note that the formation of the finger is accom- 
panied by a relaxation of the elastic strain in its vicinity, 
where the front is seen to recede (Fig. |3ja-b)). The two 
phenomena arise instantaneously - within experimental 
accuracy -, and are far more rapid than the growth rate 
of the whole pattern. In some cases, one observes small 
fingers, which start to develop and which eventually re- 
cede, if they are surrounded by larger ones (see inset (3) 
of fig. [2|center). Again this is a clear signature of the 
elastic nature of the process. 

As can be seen from the perimeter temporal evolution 
(Fig.[2]-center), the instability sets in when the strain ex- 
ceeds a critical value 7c which depends on the perimeter 
of the initial bubble Pi as 1/Pi. For a given P^, does 
not depend on the shear modulus (Fig.|4ja)), in contrast 
with the strain to fracture 7^ = (P^ — P^j/P^, which de- 
creases with G. This latter behavior tends to indicate 
that the elongational stress to failure should decrease 
with G, which is usually the case for solid brittle ma- 
terials. The characteristic wavenumber of the pattern, 
as computed from the power spectrum of the interface 
r{6) before fracture (r and 6 are the polar coordinates of 
the points lying on the pattern border) does not depend 
on G either (Fig.Qb)). 

Finally experiments performed with Setup 1, for which 
the gap b can be varied easily, allow to estimate the role 
of confinement. The width 5 of the fingers at their root, 
measured when they are just formed (Fig. [sj^a)), and the 
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FIG. 4: Dependence on the shear modulus (color online) 
(a): Critical strain at the onset of fingering, 7c, (o), and strain 
to fracture 7i^,(n), as a function of the shear modulus G 
(setup 2). (b): Power spectrum of r{0, t), averaged over time 
during stage (3) for three different bulk moduli (blue): G = 
140 Pa, (red): G = 660 Pa, (green): G = 1060 Pa (setup 2). 
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FIG. 5: Influence of the gap (color online) (a): Width of 
the fingers at their root, 5 (main panel), (b): Asymptotic 
inverse wave number n^"^(+), and ((+)) vs. b\ the 

continuous lines are indicative fits of the data, by respectively 
a power law b^^^ (main panel), b^^^ and a linear fit. 

inverse wavenumber of the instability n^^ (Fig. [sj^b)) 
exhibit nonlinear dependence on b. Indicative fits sug- 
gest power law-like behaviors 5 ~ b^^^ and n^^ b^l^ . 
Fig. |6ja) displays a spatio-temporal diagram of the local 
strain ^loc = {r{6^t) — r(6>, 0))/r(6>, 0) recorded in Setup 
1 for a gap of 2.9 mm. Quite remarkably, despite their 
sequential appearance, fingers end up regularly spaced, 
without having drifted. One can even identify the loca- 
tion of a missing finger {0 — 27r/3 in Fig. [oj^a)), which 
would eventually develop, if the pattern were not to frac- 
ture before. This suggests the existence of a well defined 
wavenumber. Fig.|6jb) shows the evolution of strain both 
along each finger, and between the fingers. The strain is 
evenly distributed until the fingers start to grow (between 
t 18s and t 25s). Then blue curves take off, showing 
stress concentration at each finger tip, while red curves 
decrease, showing that this local stress concentration is 
accompanied by a stress release in between fingers. This 
release is what causes the front to recede in the close 
vicinity of fingers, as shown above (Fig. [3|. 

Discussion We have reported experimental evidence of 
a bulk fingering instability arising in a confined layer of a 
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FIG. 6: Spatio-temporal diagram : (color online) (a): 
Local strain jioc{0,t); (b): Evolution of jioc as a function of 
time: — Blue curves: along the fingers (indicated by ticks in 
solid white lines in (a)); — Red curves: between two adjacent 
fingers, (as indicated by ticks in dotted-dashed white lines in 
(a)), b 2.9mm. 

purely elastic gel. This bulk instability is clearly different 
from the interfacial fingering observed in peeling [9l [10] 
or in some probe tak experiments For these insta- 
bilities, the wavelength increases linearly with the cell 
thickness, whereas here the inverse azimuthal wave num- 
ber ^ and the finger width (5 ~ h^l^. These 
power laws are only indicative, in view of the limited 
range of values of h which can be explored, but such non- 
linear behaviors are in clear contrast with the interfacial 
crack case. As a matter of fact preliminary probe tak 
experiments suggest that the present instability is simi- 
lar to the one reported by Shull et al. [151 [16]. Besides, 
these power laws also suggest the relevance of a second 
length scale in the problem, which should be the size I^o 
of the initial bubble. The l/I^o dependence of the critical 
strain is a strong indication in this direction. 

Given the apparent versatility of the phenomenon, one 
may wonder why such elastic fingering instabilities have 
not been characterized in details before. One key point is 
that polyacrylamide adheres to the glass plates. Blowing 
or sucking air into the cell results in considerably deform- 
ing the layers of gel which are unable to glide. Fingering 
is a way for the system to release stresses except at the 
tip of the growing finger, and thus to minimize the elas- 
tic energy. Behind the finger tip, the layer of gel is still 
pulled hard, and will eventually break at a later stage. 
The two stages - instability and fracture - may be very 
close to one another if the material is not deformable 
enough. Fracture may even occur before the instability 
can be observed. Hence the conditions for observing the 
instability is to work with adhesive, highly deformable, 
still purely elastic materials, three conditions not so often 
met, when not looked for. 

At the present stage, the precise mechanism of the in- 
stability remains unclear. However, the fact that the 
critical strain does not depend on the shear modulus sug- 
gests that there is no contribution other than elastic to 
the energy of the system. In particular, the interfacial 
tension between the injected fluid and the gel does not 



play any role. We deal with a fully three dimensional elas- 
ticity problem, with interfaces deforming on large scales, 
which requires further investigation. Moreover, the hys- 
teretic character of the phenomenon suggests that any 
linear stability analysis would be inefficient. A numeri- 
cal non-linear approach similar to the one recently used 
to understand the formation of sulci in the Biot instabil- 
ity ^ may have to be developed. As a matter of fact, 
the instability considered in [22 , which occurs when an 
incompressible elastomer is uniformly stressed laterally, 
is also purely elastic, and subcritical. The authors re- 
mark that when this instability is investigated without 
breaking translational symmetry at the surface, e.g. in 
a swollen, adhered layer of gel, then it naturally leads to 
extreme sensitivity, and subcritical instabilities connect- 
ing with buckling. And actually, the morphology of the 
patterns observed just before fracture are qualitatively 
reminiscent of buckling patterns described theoretically 
by Djondjorov and coworkers [23^. 

Finally, let us recall that fingering is commonly asso- 
ciated to liquids, while fracture is associated to solids. 
In this experiment, we show that fingering occurs also 
within the bulk of a purely elastic soft confined solid. 
Further experiments on different materials are currently 
being performed. In this context, it will be particularly 
interesting to analyze the crossover from our instability 
to a viscous fingering instability in a viscoelastic medium, 
as a function of the injected flow rate. 
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